% dat_general=[];str=[];xs=[];Disp=[];EX_general=[];DX_general=[];
% dat_ori=[];
% n=35;
% v=1;
% v=3.644551e-04
% dt=1;
% shift=200;
% 
% for i=5:10:n
%  
% Pic_Num=5000*i;
% 
% if (i==0)
%     Pic_Num=30000;
% end
% 
% 
%     dat=load(strcat('Statistical_data_concentration_X_',int2str(Pic_Num),'.sta'));
%     [nx,ny]=size(dat);
% 
% dat_ori=[dat_ori,dat];
% dat1=[];
% 
%     for j=shift+1:nx
%         dat1=[dat1;dat(j)];
%     end
%     for j=1:shift
%         dat1=[dat1;dat(j)];
%     end
%     
% dat_general=[dat_general,dat1*v*dt*Pic_Num];
% x=1:nx;
% x=x-(nx-shift);
% x=x/(dt*v*Pic_Num);
% plot(x,dat1*v*dt*Pic_Num);
% hold on 
% %dat_general=[dat_general,dat1];
% end

% x=1:nx;
% x=x-(nx-shift);
% plot(x,dat_general);


shift=200;
v=3.644551e-04;
Pic_Num=125000;
dat=load(strcat('Statistical_data_concentration_X_',int2str(Pic_Num),'.sta')); 
datss=dat;

[nx,ny]=size(dat);
sum_dat=sum(dat);
dat=dat/sum_dat;

mean_eps=v*Pic_Num;
dat1=[];

    for j=shift+1:nx
        dat1=[dat1;dat(j)];
    end
    for j=1:shift
        dat1=[dat1;dat(j)];
    end
x=1:nx;
x=x-(nx-shift);
x=x/(v*Pic_Num);
plot(x,dat1*v*Pic_Num,'x');
hold on


%=========exp bentheimer====================
a=[-0.118207, -0.004658
-0.034309, 0.042573
-0.020013, 0.08039
-0.019668, 0.113487
-0.019372, 0.141855
-0.011781, 0.203317
-0.004338, 0.250594
-0.003697, 0.312059
0.010746, 0.36406
0.03209, 0.411329
0.072758, 0.312014
0.106869, 0.250528
0.148621, 0.255231
0.225322, 0.278827
0.309122, 0.316602
0.392872, 0.34965
0.462722, 0.382705
0.539571, 0.420485
0.60242, 0.448817
0.67222, 0.477144
0.742119, 0.514928
0.825672, 0.529062
0.916324, 0.557378
1.027679, 0.571497
1.111084, 0.571447
1.215192, 0.557201
1.284549, 0.542976
1.367708, 0.519286
1.443916, 0.495601
1.513075, 0.462463
1.596086, 0.424589
1.679146, 0.391443
1.755207, 0.353573
1.851971, 0.301506
1.948833, 0.258894
2.080398, 0.211536
2.18431, 0.178378
2.27437, 0.149956
2.371479, 0.130986
2.461587, 0.107292
2.558647, 0.083594
2.655904, 0.078809
2.759914, 0.055106
2.857122, 0.045593
2.981983, 0.021878
3.148744, 0.017052
3.38496, 0.007456
];
plot (a(:,1),a(:,2),'r');
hold off
%==================================================




